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Abstract 

The challenging problems, in the field of control of chaos or of 
transition to chaos, lie in the domain of infinite-dimensional systems. 
Access to all variables being impossible in this case and the controlling 
action being limited to a few collective variables, it will not in general 
be possible to drive the whole system to the desired behaviour. A 
paradigmatic problem of this type is the control of the transition to 
turbulence in the boundary layer of fluid motion. By analysing a 
boundary layer flow for an ionized fluid near an airfoil, one concludes 
that active control of the transition amounts to the resolution of an 
generalized integro-differential eigenvalue problem. To cope with the 
required response times and phase accuracy, electromagnetic control, 
whenever possible, seems more appropriate than mechanical control 
by microactuators. 



1 Introduction 

Control of chaos or of the transition to chaos has been, in recent years, a 
very active field (see for example Ref. and references therein). Several 
techniques were developed and tested, mostly for low dimensional dynamical 
systems. The challenge lies now on finding out whether these techniques 
extend to infinite-dimensional systems. 
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A first aspect preventing a simple extrapolation of the finite-dimensional 
techniques is the fact that only a small subset of variables (or some integrated 
collective variable) is acessible to measurement. Likewise the variables on 
which one may act for controlling purposes are even more limited. A second 
aspect is that, rather than to stabilize an unstable periodic orbit (a single 
mode), what one aims in general is to suppress a continuous set of unstable 
modes, or to stabilize a particular collective mode and, at the same time, 
prevent all other modes from developing. In this sense the problem is no 
longer a standard control problem to be handled by pole placement, sliding 
mode or other standard technique. Instead, as suggested by the problem 
discussed in this paper the control problem amounts to the solution of a 
generalized integro-differential eigenvalue problem. 

A problem of both theoretical and practical importance is the control of 
the transition from laminar to turbulent motion in a boundary layer flow. 
I deal with this problem mostly as an example and prototype of the kind 
of questions and mathematical framework to be expected in the control of 
chaos for infinite-dimensional complex systems. However for the benefit of 
the reader less familiar with aerodynamical issues I have included a few 
remarks on the physical and technological context of the problem. 

By delaying the laminar to turbulent transition, an order of magnitude 
reduction in the skin friction drag is achieved. The technological benefits that 
may be derived from this reduction, led to the proposal of several methods 
for the control of the boundary layer transition. They are both of passive 
and active type and include pressure gradient control, wall suction, wall 
temperature control, polymer coating, compliant walls, etc. 

In passive type control[0 the aim is either to induce a modification of 
the curvature of the velocity profile, or to break the eddies and absorb their 
energy. 

On the other hand the active control methods, that have been proposed]^], 
aim at cancelling the growth of the Tollmien-Schlichting (TS) waves, a known 
precursor of the transition instability. This is achieved by creating a distur- 
bance of opposite phase to cancel the TS waves. The wave cancelling distur- 
bance may be created, for example, by modulated suction and blowing or by 
mechanical microactuators. This control requires an accurate set of sensors 
and actuators. The reaction time of the actuators is critical to achieve con- 
trol, especially if one aims at the feedback cancellation of nonlinear effects. 
The fact that some of the spatial growing modes have high frequencies, leads 



2 



to the suspicion that mechanical sensors and actuators, even if highly minia- 
turized, will have an hard time to deal with the high frequency instabilities 
that are known to be present in the transition. 

Greater speed and flexibility would be achieved were it possible to act on 
the flow by electromagnetic fields. With the possible exception of electrolytes 
like seawater, a direct electromagnetic action on the unmodified fluid 
im im Ell does not seem possible. However, even for neutral fluids, 



improved control of the boundary layer flow might be achieved by injecting in 
the leading edge of the airfoil a stream of ionized gas, creating a thin ionized 
layer which might then be acted upon by electromagnetic fields. In Ref. a 
detailed discussion is carried out of the effect of a streamwise directed electric 
field on the velocity profile of an ionized boundary layer, taking into consid- 
eration the fact that an injected stream of ionized gas leads to a nonuniform 
charge profile. The study establishes reference values and design estimates 
for the electric fields and ionization densities required for a significant change 
of the velocity profile. 

In the present paper a methodology is studied to assess the possibilities 
of electromagnetic control of the TS precursor waves. Usually one thinks 
of active control in terms of laminarizing the boundary layer flow. However 
the opposite situation may also occur because, for example, in stalling prone 
situations it might be useful to induce turbulence to avoid separation. Then 
the fast reaction time of electromagnetic control might also be an asset. 



2 The stability equations 

Consider the Navier-Stokes equation 

— + (f7.V)f7 = -—Vp + VAU + + X B 

Pm Pm CPm 

for an incompressible ionized fluid in an external electromagnetic field 



dt 



+ V ■ (pmU) = (2) 



In orthogonal curvilinear coordinates, denote by {u, v, w) the streamwise, the 
normal and the spanwise components of the physical velocity field U. Define 



3 



also reference quantities and adimensional variables 
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Typical values for the reference quantities, as used before[^, are Ur = 100 m 



Im, 5r = 10-^ m, pr = 1.2 Kg m-^ Er = 500 V cm'^ 



1.5 X 10"^ s-\ Then = 



15 

6.66 X 10^ and and 



/iC cm 

-j^ = 10 ~ are small quantities. 

Expressing (|lD in the adimensional variables assuming that the prod- 
uct k6 of the airfoil curvature times the boundary layer width is small and 
neglecting terms of order -R^^, and ^ one is left with 
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(4) 



(7 = 62.5 and u; = 0.15 for the reference 



where 7 = ^j^^ and u 
values above) 

The aim is to study the stability of the steady state (laminar) solutions of 
the above equations with regard to the precursor waves that develop in the 
transition region. Therefore the variables are decomposed into steady state 



[u. 



...) and fluctuating components (m , ...) 



u = u + u 

V = V + V 

w = W + w 
p = p + p 

e = e + e' 



(5) 



and one looks for normal mode solutions of the form 



u 

V 

w 

I 

p 




exp \i [ax + jSz — Vtt)} 



(6) 
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with a similar, but y— independent, form for the electric field 



E' = Eexp{i{ax + pz-nt)} (7) 

From the control point of view this implies the capability to have the elec- 
tric field react to the fluctuating velocity field with the same frequency and 
wavelength, but eventually with some delay represented by the phase of the 
complex amplitude E. To have this feedback response, a distributed set 
of sensors should be available on the surface of the airfoil. The sensors, 
of course, cannot measure the velocity field itself but only some integrated 
effect, observable at the coordinate y = (see below). 

In the transition region the quasiparallel hypothesis for the stationary 
solution is a good approximation. Namely v = w = ^ = 0. This holds for 
example for one of the scaling solutions in Ref.fl 



M = Ue 1 - exp -y—= (8) 




for X = \ and a charge distribution profile 



It is the stability and controllability of this solution that is going to be 
studied. 

By differentiating Eqs.(^ the pressure terms may be eliminated. Then, 
keeping only the linear terms in the fluctuating fields and using and (|^), 
one obtains 

vujv + %vv = la iu + u j V + ta iu + u j V 

+1^ [a [aE, + (3E,) - tf-^a {a^ + (3^) Ey] (10) 
vujw" + iOw = iau'w + iauid' — — la' Ez — iP-y-aEy} 

with the boundary conditions 

v{0) = v{oo) = w{0) = w{oo) = V (0) = V (oo) = w (oo) = (11) 
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The last three boundary conditions are obtained from the continuity equation 



iau + v + ij3w = (12) 

and the last equation in (|). 

In situations where the spanwise fluctuations may be neglected, the flow 
becomes two-dimensional and a stream function may be defined for the waves 



- dt . ' _ _d± 

oy ' ox 



(13) 



^p = (f){y) exp {i {ax + (3z — ilt)} 
Then the stability equation is 

+ i9(f)" = ia (u(f)" + u'V) - ^ (o-'-E^ - iaf^aEy^ (14) 
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which is a simplified version of the Orr-Sommerfeld equation with a driving 
term. The simplification arises from the fact that terms of order 1/Rl and 
(5^ / have been neglected. In this form the equation may be integrated once 
and reduced to a third order problem (see below). 



3 Stability and controllability results 

Consider first the spanwise stability of the scaling solution (H) without control 
{Ez = Ey = 0). The second equation in (|r^) may be integrated once and the 
integration constant set to zero using the boundary conditions ([TT|). Using 
the scaling solution (||) for u and changing coordinates to 

7] = l-exp (^-y-^^ (15) 

one obtains 

f d] 

I (1 - vf-^ ~ ^ ^^'dj]\^ ^ ^^^^ ^ iciiri''J^ (16) 

where 6i = and ai = , with boundary conditions 

w{0)=w{l) = (17) 
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Discretizing the (0, 1) interval, the calculation of the largest growing modes 
becomes an algebraic generalized eigenvalue problem which is dealt with by 
the QZ algorithm. In Fig. la one plots the largest value of Iie{iai) for real 
6i and in Fig. lb the largest value of Re(— i^i) for real ai. All modes having 
negative real parts, the conclusion is that the scaling solution is both space- 
and time- spanwise stable. Therefore we may take id = and use, for the 
streamwise stability, the stream function stability equation ([T^) . 

For the scaling solution (H), with the same change of variables, neglecting 
the term in Ey because -f^ is a small quantity, integrating once the equation 
and fixing the integration constant with the boundary conditions, the result 
is the equation 

{(1 - - 3(1 - + {l-v)i}<P + ^e,{l-r,)f^ 

with boundary conditions 

0(O) = 0'(O) = 0(1) = O (19) 

Let first = (uncontrolled equation). Using as before a finite difference 
method and the QZ algorithm with the boundary conditions imposed as three 
of the equations in the algebraic system, one obtains, for the largest value 
of Re(2ai) and real 6i, the results shown in Fig.2. It means that there is 
a range of frequencies for which there is spatial growth of the streamwise 
fluctuations. Therefore the scaling solution is spatially unstable. 

To derive a controlled equation one has to realize that the only physical 
quantities, that it is reasonable to assume to be observable, are the pressure 
fluctuations on the airfoil or the integrated effect of the electrical current as 
seen at the surface of the airfoil. Pressure fluctuations may be detected by a 
distributed set of microphones and the integrated electrical current fluctua- 
tions are essentially the induced magnetic field fluctuations on the spanwise 
direction. For definiteness I will assume that a set of sensors is available to 
measure the effect of the electrical current fluctuations. To achieve control 
this measurement is used to modulate a variable component of the applied 
electric field. That is 

E, = k dyao (l - ^) u{y) (20) 
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with k a complex proportionality constant, the meaning of the phase being 
the control delay. Then, the controlled equation is 

{(1 - v)% - 3(1 - v)% +{i-v)ii}<f>+Mi-vri 

with C = ^^^^ . The controllability problem amounts to find out whether 
all eigenvalues have negative real parts in the integro-differential problem 
defined by Eq.(pll). Let 6^1 = 60 , the frequency for which the largest Re(zai) 
is at its maximum. For real C, Fig.3 shows that for C > 1 the largest mode 
has spatial decay, hence the solution becomes stable. For the results in Fig. 4, 
let C = \C\ e**^ with \C\ = 1.5 and variable phase if. One sees that there is 
a range of phase delays which stabilize the solution and, conversely, outside 
this range the solution is strongly unstable. 



4 Conclusions 

1. Stabilization of a stationary configuration in an infinite-dimensional sys- 
tem involves the study of infinitely many disturbance modes, some of which 
may grow in time and space. In addition, the measurable observables, to 
which some local control may react, involve the integrated effect of many 
variables. Therefore the mathematical structure of the problem to be solved 
is expected to be, as in this example, an integro-differential generalized eigen- 
value problem. 

2. The unstabilizing disturbances that need to be controlled in extended 
systems have in general a nontrivial space-time structure. Therefore a set of 
distributed sensors and actuators is needed to achieve a space-time control- 
ling action. 

3. The fact that, in practice, only global integrated variables are ob- 
servable, restricts the feedback control to these variables only. Therefore, 
for extended systems, there is no guarantee that control will be achieved in 
general and success is only to be expected in particular favorable cases. 

4. The laminar to turbulent transition, in the boundary layer, begins 
with the appearance of downstream moving waves which at first grow slowly 
and may be described by a linearized equation. After reaching a certain 
amplitude however, the waves develop strong three-dimensional structures 
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and nonlinearities and a rapid transition to turbulence becomes unavoidable. 
Therefore, if effective control is to be achieved, it is essential to have a fast 
and locally accurate feedback to tame the instabilities before the spanwise 
differential amplification of the TS waves begins to occur. Because it is prob- 
ably very difficult to obtain the required speed and accuracy with mechanical 
microactuators, electromagnetic controlling schemes seem worth to explore. 

5 Figure captions 

Fig.l - (a) Space stability of the spanwise modes, (b) Time stability of the 

spanwise modes 

Fig. 2 - Space instability of the streamwise scaling solution 

Fig. 3 - Largest Re(zai) for the controlled equation (6*1 = 60 and real C) 

Fig. 4 - Largest Re(zai) for the controlled equation {6i = 60, \C\ = 1.5 

and variable phase) 

References 

[1] T. Shinbrot; "Progress in the control of chaos". Advances in Physics 44 
(1995) 73. 

[2] D. M. Bushnell and J. N. Hefner (Eds.); "Viscous drag reduction in 
boundary layers". Progress in Astronautics and Aeronautics, vol. 123, 
American Institute of Aeronautics and Astronautics, Washington 1990 
and references therein. 

[3] R. Vilela Mendes and J. A. Dente; "Boundary-layer control by electric 
fields: A feasibility study" , Report no. [physics/ 9705020| . 

[4] R. D. Joslin, G. Erlebacher and M. Y. Hussaini; "Active control of insta- 
bilities in laminar boundary layers - Overview and concept validation" , 
Journal of Fluids Engineering 118 (1996) 494 and references therein. 

[5] H. Schlichting; "Boundary-Layer Theory" 6th. edition, MacGraw Hill, 
New York 1968. 

[6] A. D. Young; "Boundary Layers", BSP Professional Books, Blackwell, 
Oxford 1989. 



9 



[7] A. K. Gailitis and O. A. Lielausis; "On the possibility of drag reduction 
of a flat plate in an electrolyte", Appl. Magnetohydrodynamics, Trudy 
Inst. Fis. AN Latv. SSR 12 (1961) 143. 

[8] A. B. Tsinober and A. G. Shtern; "Possibility of increasing the flow 
stability in a boundary layer by means of crossed electric and magnetic 
fields", Magnetohydrodynamics 3 (1967) 103. 

[9] H. K. Moffat; "On the suppression of turbulence by a uniform magnetic 
field", J. Fluid Mech. 28 (1967) 571. 

[10] C. B. Reed and P. S. Lykoudis; "The effect of a transverse magnetic 
field on shear turbulence", J. Fluid Mech. 89 (1978) 147. 

[11] A. Tsinober; "MHD flow drag reduction", in 0, page 327. 

[12] C. Henoch and J. Stace; "Experimental investigation of a salt water 
turbulent boundary layer modified by an applied streamwise magneto- 
hydrodynamic body force", Phys. Fluids 7 (1995) 1371. 

[13] C. H. Crawford and G. E. Karniadakis; "Reynolds stress analysis of 
EMHD-controUed wall turbulence. Part I. Streamwise forcing", Phys. 
Fluids 9 (1997) 788. 



10 



Fig. 1a 
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